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A recent pubhcation [Phys. Rev. Lett. 100, 140408 (2008)] shows that there is a relation between 
mean-field chaos and multi-particle entanglement for BECs in a periodically shaken double well. 
"Schrodinger-cat" like mesoscopic superpositions in phase-space occur for conditions for which the 
system displays mean-field chaos. In the present manuscript, more general highly-entangled states 
are investigated. Mean-field chaos accelerates the emergence of multi-particle entanglement; the 
boundaries of stable regions are particularly suited for entanglement generation. 
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I. INTRODUCTION 



results for the emergence of entanglement. 



Systems of ultra-cold atoms with time-periodic poten- 
tial differences have recently started to be implemented 
experimentally both on the single-particle level PJ and 
for Bose-Einstein condensates (BECs) [2] [3]. These ex- 
periments realise schemes of tunnelling-control devised 
in recent years [U [S] . Very recently, even the control of 
the Mott-insulator transition via periodic shaking pi 171 
has been realised experimentally 8 . This would, roughly 
speaking, correspond to making a glass of water freeze by 
turning on a loud-speaker. Proposals for shaking-induced 
tunnelling control for ultra-cold atoms not realised so far 
include schemes for entanglement generation [Tni ITT] . 
The tunnelling control based on shaking the potentials 
uses concepts successfully implemented in the dressed- 
atom picture [T^. 

The double-well potentials [ISj used for the single- 
particle experiments [1] would promise further excit- 
ing experimental results for BECs: a simple mean-field 
model of this system is known to be chaotic for not too 
small interactions [H [TH [TSJ [HI [TS]. For BECs in a 
double-well potential with periodic shaking, we demon- 
strated numerically the emergence of "Schrodinger-cat" 
like states in phase-space if the parameters are chosen 
such that the mean-field system displays a coexistence of 
regular and chaotic dynamics |19j . In this paper, more 
general entangled states are investigated. As shown by 
the authors of Ref. [20 , the quantum Fisher informa- 
tion j5TJ [22] can be used to identify multi-particle en- 
tanglement. This leads to a a greater variety of highly 
entangled states than we investigated in Ref. TOJ. 

This manuscript is organised as follows: Section |IT] in- 
troduces the models used to describe the BEC in a double 
well both on the mean-field (Gross-Pitaevskii) level and 
on the TV-particle quantum level. Section [nT| explains the 
way entangled states are detected in the numerics. Be- 



fore concluding the paper. Sec. IV shows the numerical 



II. MODEL 

Bose-Einstein condensates in double-well potentials 
are often described via a model originally developed in 
nuclear physics |23j , a multi-particle Hamiltonian in two- 
mode approximation |24j : 



H 



-I- /i(^o + Ml sin(wt)) ^0^02 — ojo^^ , (1) 



where the operator o^-^-* creates (annihilates) a boson in 
well j; hfiQ is the tilt between well 1 and well 2 and hfii 
is the driving amplitude. The interaction between a pair 
of particles in the same well is denoted by 2hK. Appli- 
cations of such Hamiltonians include multi-particle en- 
tanglement [TTl , high precision measurements, 
many-body quantum coherence \27\ I28j and spin sys- 
tems P] . 

Within the approximation on which the Gross- 
Pitaevskii equation is based, a multi-particle wave- 
function is characterised by two variables, 6 and (f). The 
fraction of atoms found in well 1 (well 2) is given by 
cos^ [6 1 2] (sin'^ [^/2] ) ; the phase- factor between both wells 
is given by exp(i<^). All atoms occupy this single par- 
ticle wave- function. One can also find corresponding 
TV-particle wave-functions known as "atomic coherent 
states" [3D]. In an expansion in the Fock-basis |n, TV — n) 
with n atoms in well 1 and N — n atoms in well 2 these 
wave-functions read: 



^ /TV^ 
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The Gross-Pitaevskii dynamics can be mapped to that 
of a nonrigid pendulum [31] . Including the term describ- 
ing the periodic shaking, the Hamilton function is given 



2 



used to define 



N 






0.5 


Jj. — 


1 f.u 

- [a{a^ + a'^a^ 


1 




Jy = 









■h = 





Q 



FIG. 1: (Colour online) An ideal "Schrodinger-cat" like 
state which is the superposition of two atomic coherent 
states (|2| with hardly any overlap. The figure shows the 
projection \{'ipca.t\d,4>)\'^ of the "Schrodinger-cat" like state 
V'cat = -0.6, (j> = 1.2) + \z^ 0.65, <t) = -2.20)) onto 

the atomic coherent states \6,(t>) (e.q. [2| in dependece of 
z and (j). While Ref. [TS] concentrated on identifying such 
highly entangled state, the present manuscript uses a differ- 
ent approach to search highly entangled states: the sufficient 
condition (|9| derived in Ref. by using the quantum Fisher 
information Q (Ref. [20] and references therein). 



by (z = cos^{9/2) - sm^{0/2)): 



Nk 
2z 



\J\ — cos((/)) 



= m 



(3) 



In our case, 2;/2 is the experimentally measurable |13j 
population imbalance which can be used to characterise 
the mean-field dynamics. The dynamics for this system 
are known to become chaotic [32] . 



III. MULTI-PARTICLE-ENTANGLEMENT & 
QUANTUM FISHER INFORMATION 

Multi-particle-entanglement ^51 1^ ^5] |M1 [55] [55] [571 
[38] is a hot topic of current research; to experimentally 

realise "Schrodinger-cat" like superpositions of BECs is 
still a challenge of fundamental research. For the period- 
ically driven double-well potential, the relation between 
emergence of "Schrodinger-cat" like mesoscopic super- 
positions in phase-space and mean-field chaos has been 
discovered in Ref. [13] • An ideal example of such a meso- 
scopic superposition is shown in Fig. [T| the fidelity [S3] 
of some of the highly entangled states found numerically 
was well above 50% |19j . In this manuscript, we employ 
the fact that the quantum Fisher information can be used 
to detect multi-particle entanglement [2U] . 

Before defining the quantum Fisher information, we 
note that the creation and annihilation operators can be 



(4) 
(5) 
(6) 



which satisfy angular momentum commutation rules. 
Except for a factor of A^, the operator is the oper- 
ator of the (experimentally measurable [131 ) population 
imbalance. 

While the quantum Fisher information i^q can be de- 
fined for statistical mixtures, it is particularly simple for 
pure states [20] : 



Fq = 4(A J„)2 ; 



(7) 



where A J„ are the mean-square fiuctuations of J in di- 
rection n. In the following we choose the z-direction, 
thus 



(8) 



For iV particles, a sufficient condition for multi-particle 
entanglement is given by [20]: 



Font > 1 



N 



For the ideal "Schrodinger-cat" state in real space, 
1 



noon; 



V2 



(|7V,0) + |0,iV)) 



(9) 
(10) 

(11) 



this entanglement fiag reaches a value of N. Thus, while 
Eq. ^ already indicates multi-particle entanglement, 
values of 

Font » 1 (12) 

demonstrate highly entangled states. 



IV. RESULTS 

In order to characterise whether or not the mean-field 
dynamics are chaotic, Poincare surface of sections are 
particularly suited: for a set of initial conditions, the 
mean-field dynamics of the periodically driven system 
characterised by the angular frequency lu is plotted at 
integer multiples of 27r /lu. Figure 2] shows a phenomenon 
which is very characteristic for the non-rigid pendulum 
on which the BEC in a double well was mapped within 
mean-field: the coexistence between chaotic and regular 
regions. 

While each initial condition will, in general, lead to 
many dots in plots like Fig. [2] we proceed in the spirit 






FIG. 2: Poincare surface of section for the mean-field sys- 
tem. Closed loops characterise stable orbits whereas chaos 
is represented by irregular dots. The parameters are chosen 
such that they correspond to a one-photon resonance: a tilt 
of 2fio/^ = 3.0, a driving frequency of a; = 3Q, an interaction 
of Nn/Q. = 0.8 and a driving amplitude of 2^i/fJ = 0.9 (cf. 
Ref. 0). 




of Ref. [T5] to characterise the A^-particle dynamics: In 
Fig. [3] each initial condition leads to only one point: the 
entanglement flag (10 1 for this initial condition after a 
fixed time tVl. Highly entangled states (Eq. (12l) can 



be found on the boundary of stable regions; already for 
short times (Fig. [s] a) such highly entangled states can 
occur. For larger times (Fig. [s] b) many features of the 
Poincare surface of section in Fig. [2] are visible in the 
entanglement generation which essentially spreads over 
the entire chaotic part of the initial conditions. 

For parameters which display no chaotic parts in the 
Poincare surface of section (Fig. |4|, on short time- 
scales hardly any entanglement emerges (Fig. |5]a). For 
larger time-scales, entanglement generation occurs on the 
boundaries of stable regions (Fig.|5]b). 



V. CONCLUSION 

In Ref. [12] we discovered the relation between meso- 
scopic "Schrodinger-cat" like superpositions in phase 
space and mean-field entanglement. In the present pa- 
per, we demonstrated that also for more general entan- 
gled states, multi-particle entanglement can be a quan- 
tum signature of chaos. For regular systems, the general 
entangled states also occur. However, it is restricted to 
the boundaries of stable regions and only occurs on longer 
time-scales. While the focus in Ref. [12] was on finding 
particularly highly entangled states for each initial con- 
dition, in the present manuscript the entanglement pro- 
duction was shown for each initial condition at the same 
point in time, both for short and longer times between 



FIG. 3: (Colour online) The entanglement flag in a two- 
dimensional projection as a function of the initial condition 
for A'^ = 100 particles. All other parameter as in Fig. [2] The 
TV-particle wave-function which corresponds to the mean-field 
initial conditions (zo,<^o) can be found in Eq. ([2|. In the 
upper panel, the (experimentally measurable) entanglement 
flag (10 1 is shown after a time of — 10, in the lower panel 
the time is tfl = 100. Highly entangled states (cf. Eq. (12l) 
occur in the entire chaotic regime. 
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FIG. 4: Poincare surface of section for the mean-field sys- 
tem (cf. Fig. [2|. The parameters are chosen such that they 
correspond to a 3/2-photon resonance [5! with Nn/Q. = 0.1, 
2^0/n = 3.0, u/Q. = 2.08 and 2^i/Q. ='l8. 
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the onset of the computer experiment and the read-out. 



As an entanglement-flag, we apply the quantum Fisher 
information for pure states. In this paper we use it in a 
way which is particularly easy to measure experimentally. 



However, using Eq. (10 1 assumes a pure state and would 
thus only be valid in an ideal system without decoher- 
ence. For more realistic situations, experimental signa- 
tures as in Refs. [501 SI] should be investigated in the 
future. 




FIG. 5: (Colour online) The entanglement flag ( 10 1 in a 



two-dimensional projection as a function of the initial con- 
dition (cf. Fig. Jsl for A'^ = 100 particles; all other param- 
eter as in Fig. Hf For the regular regime given by the pa- 
rameters of Fig. |4j after short time-scales of tQ, = 10 (upper 
panel) entanglement-production is very weak (note that the 
brightness-entanglement-coding difi'ers by a factor of 10 from 
Fig. [5] and the longer time-scales til = 100 depicted in the 
lower panel). For larger time-scales, entanglement produc- 
tion mainly occurs on the boundary of stable regions (lower 
panel). 
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